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The measurements of the isoscalar giant monopole resonance (GMR), also called the breathing 
mode, are analyzed with respect to their constraints on the quantity Mc, e.g. the density dependence 
of the nuclear incompressibility around the so-called crossing density pc—Q-l fm~^. The correlation 
between the centroid of the GMR, Eqmb., and M^ is shown to be more accurate than the one 
between Bgmr and the incompressibility modulus at saturation density, K^o, giving rise to an 
improved determination on the nuclear equation of state. The relationship between M^ and Ka^ 
is given as a function of the skewness parameter Q^o associated to the density dependence of the 
equation of state. The large variation of Qoo among different energy density functionnals directly 
impacts the knowledge of K^o'- a better knowledge of Qoo is required to deduce more accurately 
Koo ■ Using the Local Density Approximation, a simple and accurate expression relating -Egmr and 
the quantity Mc is derived and successfully compared to the fully microscopic predictions. 

PACS numbers: 21. 10. Re, 21.65.-f, 21.60.Jz 



I. INTRODUCTION 

The determination of the nuclear incompressibility is 
a long standing problem. The earliest microscopic anal- 
ysis came to a value of ^00=210 MeV [l[, but with the 
advent of microscopic relativistic approaches, a value of 
iiroo=260 MeV was obtained [2|- The fact that K^o can- 
not be better determined than 230± 40 MeV, taking into 
account the whole data on the isoscalar Giant Monopole 
Resonance (GMR) as well as the various methods to re- 
late the GMR to Kao (see e.g. Ill-Q) lead to a recent 
effort to reanalyse the method 'To|- 

Pairing effects and similarly the shell structure ef- 
fects on the nuclear incompressibility were analyzed along 
these lines. Since the first investigation [11[, several stud- 
ies have shown that pairing effects have an impact on the 
determination of Koo [7|, |8| , and it was considered as a 
possible cause of the difficulty to accurately constrain 
Koa ■ This effect of pairing on the incompressibility mod- 
ulus has also been analyzed in nuclear matter, showing 
that the main effect is occurring at subsaturation densi- 
ties [141 . However there is a general consensus between 
the various microscopic models that pairing effects on 
Koo are not strong enough to explain the lack of accuracy 
in the determination of the nuclear incompressibility [3- 
0, [lj|. Other effects have to be investigated. 

Recently, the density dependence of the nuclear incom- 
pressibility was re-investigated suggesting that the corre- 
lation between the centroid of the GMR and the incom- 
pressibility modulus Koo at saturation density is blurred 
by the density dependence of the nuclear equation of state 
in different models 10]. The observed differences in the 
extraction of K^o from the i^cMR are based on different 
models and attributed to the density dependence of the 
equation of state (EoS) which has still to be better con- 
strained. The observation of a crossing point provided 
a possible path to be investigated. The crossing point 



arises from Energy Density Functionnals (EDFs) that 
are designed to describe finite-nuclei observables: their 
density-dependent incompressibility K[p) crosses around 
the mean density in nuclei, pc — 0.1 fm~'^. It was there- 
fore proposed that the quantity Mj,, e.g. the density 
dependence of K{p) around the crossing density pc is the 
quantity that shall be constrained by measurements of 
the GMR, instead of K^o- 

The aim of the present article is to further analyze 
the correlation method on which is extracted the incom- 
pressibility modulus Koo, and to give a better basis on 
the alternative method based on the correlation between 
Eqmk and Mc- A comparison between the two methods 
is given in Sec. [TTl for ^ospj^ g^^^ ^^^Sn nuclei, showing 
the relevance of the new method [10] ■ In Sec. IIIII the 
source of uncertainties in the determination of Koo is di- 
rectly related to the skewness parameter Qoo- The skew- 
ness parameter gives indeed the present limitation on the 
knowledge of the density dependence of the nuclear EoS 
between the crossing and the saturation densities. The 
origin of the crossing density is also demonstrated in the 
case of the Skyrmc EDFs. In Sec. lIVi the explicit relation 
between the centroid of the GMR and the quantity M^ 
is derived using the Local Density approximation (LDA), 
and keeping as much as possible analytical relations be- 
tween the various quantities, in order to facilitate their 
interpretation. The results are compared to the fully mi- 
croscopic one. Conclusions are given in Sec. |V| 



II. THE MICROSCOPIC APPROACH 

In this section, we first summarize the Constrained 
Hartree-Fock-Bogoliubov approach (CHFB) used to ac- 
curately predict the Isoscalar Giant Monopole Resonance 
(GMR) energy. We then provide the definition of the 
parameter Mc, driving the density dependence of the in- 



compressibility around the crossing point. Finally, using 
these two quantities, the correlation analysis between the 
GMR energy and the incompressibility modulus Koa on 
one hand, and the GMR energy and the parameter Mc 
at the crossing density on the other hand, are compared. 



A. Microscopic calculation of the GMR energy 



CHFB solution |A), 



m_i = -- 



^(AigiA) 



(6) 



A=0 



B. Constraints on the equation of state deduced 
from £gmh 



We first recall how to predict the Eqmr- We use the 
sum rule approach in order to microscopically calculate 
the centroid energy of the GMR. In such a microscopic 
approach, we calculate the energy as 



E, 



GMR — 



mi 



(1) 



where the fc-th energy weighted sum rule is defined as 

mfc=^(^,)'l(*IQ|0)P, (2) 

i 

with the RPA excitation energy Ei and the isoscalar 
monopole transition operator. 



Q = E^ 



(3) 



The calculations using fully microscopic approaches 
based on EDF are usually performed using CHFB or the 
RPA approach [ij]. In the present case we calculate the 
GMR centroid for the Skyrme EDF with the GHFB ap- 
proach. For completeness, results using other function- 
nals such as Gogny and relativistic functionnals will also 
be given. The CHFB method is known to provide an 
accurate prediction of the GMR centroid. 

In the following the energy weighted moment mi and 
the r7i_i moment are directly evaluated from the ground- 
state obtained from Skyrme CHFB calculations. The 
moment mi is evaluated by the double commutator using 
the Thouless theorem 1 1511: 



2h^A 



nil = 



ir')- 



(4) 



where A is the number of nucleons, m the nucleon mass 
and (r^) is the rms radius evaluated on the ground-state 
density given by Skyrme HFB. 

Concerning the evaluation of the moment m,_i, the 
constrained-HFB approach is used. It should be noted 
that the extension of the constrained HF method [3, 16 1 
to the CHFB case has been demonstrated in Ref. 17 1 
and employed in Q. The CHFB Hamiltonian is built 
by adding the constraint associated with the isoscalar 
monopole operator, namely 



Hr, 



^H + XQ, 



(5) 



and the m_i moment is obtained from the derivative of 
the expectation value of the monopole operator on the 



Next, the parameter Mc is defined. Instead of corre- 
lating Eqmr and Koo, it was proposed that the energy 
of the GMR ([T]) gives a strong constraint on the quan- 
tity Mc defined, at the crossing density pc — 0.1 fm~^, 
as [13, 



Mc = 3pcK'{p)\p-_ 



(7) 



where the density-dependent incompressibility K{p) is 
derived from the thermodynamical compressibility x(p) 



as m 



Kip) 



9p 



18 



x(p) 

and the pressure is 



Pip) + 9p- 



d^Eip)/A 
dp^ 



P(.)..^^ 



(8) 



(9) 



The parameter Mc was introduced instead of Kao = 
Kipo) (where po is the saturation density) in the cor- 
relation analysis based on Eqmr since i) the crossing 
density pc is closer to the average density in finite nuclei 
than the saturation density po , and ii) the crossing of the 
incompressibility at pc makes the -Bgmr mostly sensitive 
to the derivative of the incompressibility at the crossing 
density [l3| • It should be noted that the existence of a 
crossing density for other EoS quantities, such as for in- 
stance the symmetry energy [131, ^^^ neutron EoS i2(]ll2l| 
and the pairing gap in nuclear matter |22:] . was also ob- 
served. It might reveal the general trend that the experi- 
mental constraints drive these quantities towards a cross- 
ing point at around the average density of finite nuclei. 
Various EDF's shall however exhibit various density de- 
pendencies around the crossing point. At first order the 
derivative of the incompressibility (or symmetry energy 
or pairing gap) at this point will differ between various 
EDF's and additional measurements in nuclei shall char- 
acterize these derivatives. For instance, the derivative of 
the neutron EoS around pc — 0.11 fm"'^ was found to be 
strongly correlated to the neutron skin in ^"^Pb |20l.[2]|. 
giving a strong support to improved experimental mea- 
surements of this quantity 23 1 . 

Fig. [T] depicts Kip), between half of the saturation 
density and the saturation density, for several Gogny, 
Skyrme and relativistic EDFs. A large dispersion is 
observed at saturation density (p/po — 1) whereas at 
p/ Po — 0.71 there is a much more focused area, defining 
the crossing density pc- 
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TABLE I. Evaluation of Kc and of the two terms defining 
tlie incompressibility K(p) (Eq. (|8])), at the crossing density 
pc — 0.71po and at the saturation density po for a set of 
different Skyrme EDFs. 



FIG. 1. (Color online) EoS incompressibility K(p) calcu- 
lated with various relativistic and non-relativistic function- 
nals. The Skyrme EDF are in solid lines 



We now analyze the two contributions to the density- 
dependent incompressibility, as depicted by Eq. ([S]). The 
first term of the r.h.s in Eq. ^ is proportional to the 
pressure P{p), which is indeed related to the first deriva- 
tive of E/A and the second term is the second derivative 
of the binding energy E/A with respect to the density. 
The former vanishes at saturation density, by definition. 
Table U displays the expectation values of these two con- 
tributions to the incompressibility K{p) (Eq. ([8|)) at the 
crossing density pc ~ 0.71po and at the saturation den- 
sity for several EDFs. The total value of K{pc) = Kc 
at the crossing density is also displayed, while at po, 
Koo — 9/9o — g 2 I pp. The very weak dispersion of Kc 
as a function of the EDFs is striking, whereas the incom- 
pressibility modulus at the saturation density K^o is more 
scattered. At the crossing density, the incompressibility 
Kc is given as the sum of the first and second derivatives 
of the energy per particle E/A, which act in opposite 
signs, see Table ID The contribution of the pressure at pc 
is not negligible, at variance with its contribution at p^, 
and largely contributes to the stabilisation of Kc- The 
correlations between -Egmr and the solely second deriva- 
tive of E/A at the saturation density might not be the 
most appropriate one and the EDF-invariant property of 
the crossing point {pc,Kc) shall be useful. 



C. {Egmi\.,Koo) versus {Egmb.,Mc) correlation 
analysis 

Using Eqmk and Mc discussed in the previous sections, 
it is possible to determine if Mc is better constrained by 
Eqmr than Koo- The correlation diagrams (-E'cMRj-f'^cx)) 
and [Eqmk^Mc) are compared on Fig. [5] in the case of 
208p|2_ The relativistic DDME2 interaction in the cor- 



relation graph (-BgmRj^oo) is largely deviating from the 
others, as it is well known [2tJ] while it is much more 
compatible with the others in the {Eqm^^Mc) correlation 
graph ^IQj . On the contrary, restricting to the Skyrme in- 
teractions, the quantities (i?GMR,^oo) and (ii^GMR,Afc) 
are equally well correlated. This is directly related to 
the good correlation between {Mc,Koo) due to a simi- 
lar density dependence among the Skyrme EDFs (in p"), 
which will be discussed in section IILC. However, consid- 
ering various models with different density dependencies, 
a better correlation is observed between -Egmr and A'Ic^ 
compared to the one between Sgmr and K^o- It should 
be noted that in the (i?GMR,Afc) correlation graph, an 
ordering between Skyrme and relativistic models is also 
observed. 

Fig. [3] displays the same correlations in the case of the 
^^'^Sn nucleus. In this case pairing effects are known to 
slightly impact the position of the GMR |13| , leading to 
a larger dispersion compared to the ^°^Pb case. How- 
ever similar conclusions can be drawn, namely i) Mc is 
a better correlated quantity with Sgmr than K^o, ii) a 
good (SgmRj^oo) correlation is also observed among the 
Skyrme interactions and iii) the (i?GMR,-Wc) correlation 
exhibits an ordered correlation between the Skyrme and 
the relativistic EDFs. 

To conclude, these results on Mc provide a step to- 
wards compatible results between Skyrme, Gogny and 
relativistic approaches [10]. The extracted value for the 
quantity Mc in ^^°Sn and 208p|-, n^dej ^re also in bet- 
ter agreement between each other than the corresponding 
Koo values: considering the various EDFs as well as the 
i20Sn and the ^"Spb data, one gets Mc=1100 ± 70 MeV 
(6% uncertainty), and Xoo=230 ± 40 MeV (17% uncer- 
tainty) [l3| • It should be noted that these considerations 
on the slope of the incompressibility Mc at the crossing 
point have recently been used in Ref. [2^ where a good 
linear correlation between Sgmr and Mc is also found, 
including the so-called BCPM functionnal. 

In summary, using microscopic approaches, it is ob- 
served that the correlation between Mr and the centroid 
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FIG. 2. Centroid of the GMR in 2°**Pb calculated with the 
microscopic method (see text) vs. the value of Koc (top) and 
Mc (bottom) for various functionals 0, i, 53, |24-^]. The 
solid and dashed lines correspond to fit on the Skyrme and 
relativistic EDF values, respectively. 



Eqmr is less dispersive, and therefore more universal 
among various models, than the one between K^o and 
^-GMR [10|. In the next section, we shall provide a more 
quantitative understanding of the differences between the 
quantities Mc and Ko^, explaining the role of the density 
dependence of the equation of state between the crossing 
and the saturation densities. 



III. DENSITY EXPANSION OF THE 
EQUATION OF STATE 

The striking stability of Kc among the various Skyrme 
EDFs (Table H]) deserves an investigation. In this sec- 
tion, the density dependence of the equation of state is 
discussed in terms of the derivatives of the EoS with re- 
spect to the density. The stability of Kc as well as the 
relation between the slope of the incompressibility mod- 
ulus Mc and the parameters Koa and Qoo are derived, 
providing an explanation for the difficulty to constrain 

Kr^. 



A. Density dependence of the equation of state 
around po 

We start from a systematic expansion around the sat- 
uration density po of the binding energy, such as in the 
Generalized Liquid Drop Model (GLDM) [H [U where, 
in symmetric matter, the energy per particle reads 



FIG. 3. Same as Fig. [2]for i2°Sn. 



with X = {p — po)/{3po), Po being the saturation density 
of symmetric nuclear matter. Qoo is the third derivative 
of the energy per particle. 

Applying Eqs. ^ and (O to the expansion Eq. p^ . 
one obtains the pressure. 



P(x) = i(l + 3a;)2 



-"-oo*^ ~r -tc/oo-^ 



(11) 



and the incompressibility, 

K{x) ^{1 + 3x) [Xoc + (Sifoo + Qoo)x + 6Qooa;^ 



1 



E{x) ^Eao + T^KooX^ 



(10) 



(12) 



Fig. |4] displays the binding energy Eq. ((T0| , pressure 
Eq. ([Tl]) and incompressibility Eq. (|12p as function of 
the density p going from to 0.2 fm~^, for typical values 
for the quantities : -Boo=-16 MeV, Kao^2A{) MeV and 
Qoo=-350 MeV. The different curves correspond to vari- 
ous approximations in the density expansion of the bind- 
ing energy Eq. (ITUl) . For instance, the solid line in the 
binding energy E/A corresponds to the 0-th order in the 
density expansion Eq. (jlOp where only the quantity -Boo 
has been included, all other quantities being set to zero. 
The dotted-line (E-l-K) takes into account the quantities 
Eao and K^o, and the dashed line (E-l-K-l-Q) includes 
the quantities Eoo, Kqo and Qoo- Similar approximations 
have been performed in the case of the pressure Eq. (jlip 
and incompressibility Eq. (J12p . A good convergence is 
found when successively including in the expressions for 
the binding energy, the pressure and the incompressibil- 
ity the quantities Soo, Koo and Qoo- These quantities 
describe the density dependence of the equation of state 
and are given in table |ll] for a set of models considered 
in this work. 

It is clear from Table |TT] that while the quantities E^o 
and Koo are not varying by more than 20%, the values 
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FIG. 4. (color online) a) Binding energy EjA in MeV, b) 
pressure in MeV frti"^ and c) incompressibility K in MeV, as a 
function of the density for various truncation in the expansion 
Eq. (|10|l . See text for more details. 
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fm-3 


MeV 
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TABLE II. Parameters appearing in the density expansion of 
the binding energy E/A Eq. (|10p for a set of models con- 
sidered in this work: po is the saturation density, Eoo the 
binding energy, Koo the incompressibility modulus, and Qoo 
the skewness parameter. 



for the skewness parameter Qoo is almost unconstrained 
and can vary by more than 100% among the models. The 
uncertainty in the determination of the skewness param- 
eter Qoo gives, in a quantitative way, the main lack of 
knowledge in the density dependence of the equation of 
state. The uncertainty on Qoo also impacts the density 
dependence of the pressure and more, interestingly here, 
of the incompressibility. 
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FIG. 5. Koo versus Qoo for several models. The solid hne 
corresponds to the fit on the Skyrme EDFs values. 



B. Stability of Kc 

Let us now provide an explanation for the stability of 
Kc observed in Table HIl From Eq. (IT^ . and assuming the 
validity of a density expansion from po to pc, we obtain 



/C 



(1 + 3Xc) [(1 + 9Xc)Koo + (1 + 6Xc)Xc 



(13) 



with Xc = {pc - po)/(3po)- 

In the case of Skyrme interaction, there is a good cor- 
relation among the quantities Koo and Qoo, as shown 
in Fig. [SJ The parameters K^^ and Qoo are mostly de- 
termined by the same term, the term ts in p", in the 
case of Skyrme interaction ]1(|. Due to their similar 
density dependence (in p"), the Skyrme EDFs exhibit 
indeed a linear correlation among these two quantities 
whereas the picture is blurred when considering at the 
same time Skyrme and relativistic EDFs. The linear cor- 
relation among the Skyrme EDFs can be described by, 



Koo = a + bQo 



(14) 



with a = 338 ±9 MeV and b = 0.29 ±0.03. Injecting H 
in Eq. IT^ . one gets 



Kc 



{l + 3xc)[{l + 9xc)a + f{xc)Qc 



(15) 



with f{x) = (6a;2 + (96 + l)x + b). An EDF-almost- 
independent value of K^. is therefore obtained for f(x)=0 
since Qoo is the only EDF-dependent quantity in Eq. 
p^ : the solution of f(x)=0 shall therefore provide the 
crossing point observed on Fig. [TJ The function f{x) 
has only one zero for positive densities, given by Xc — 
-0.095 ± 0.002, which corresponds to p = (0.714 ± 
0.005)po = Pc- 

In the Skyrme case there is therefore a density, pc, for 
which the incompressibility modulus K{pc) is indepen- 
dent from the quantities Koo and Qoo defined in Eq. (jlOp . 



and is 



(l + 3xc)(l + 9a;c)a: 



Po Po 



-2)a. (16) 



Taking the value for a given by the hnear correlation, one 
finds Kc — 34 ± 4 MeV, confirming the value of the cross- 
ing point shown in Table |T] and on Fig. [T]for the Skyrme 
EDFs. This approach confirms in a both quantitative 
and qualitative way the existence of a crossing point, es- 
pecially in the case of the Skyrme EDF. In the case of 
the other EDFs, it is rather a crossing area that is ob- 
tained (Fig. [ij, due to their various density dependence, 
as discussed above. 



C. Relation between Kaa and Mc 

Fig. [6] displays the {Mc,Koo) correlation for four 
Skyrme EDFs. Adding to this correlation graph EDFs 
with other density dependences, such as the relativistic 
one, drastically blurs this correlation. This shall orig- 
inates from the large uncertainty on the value of the 
skewness parameter Qoo among the EDFs discussed in 
the previous section (Table Hi]). Using Eqs. ([71), dH) and 
P^ . the quantity Mc can be expressed as, 



Mr. 



iKc + (1 + 'ixcf 9Koo + (1 + 12a;c)Qoo ■ (17) 



The correlation between M^ and K^o depends on the den- 
sity dependence of the binding energy reflected in the 
skewness parameter Qoo, which can vary to a large ex- 
tent, see Table HIl More precisely, from Eq. ([TT]). one can 
deduce the value of the quantity Koc as. 



K^ 



1 Mr - 3Kr 1 + 12xc 



9(H-3a;c)^ 



9 



(18) 



where the second term of the r.h.s shows the theoretical 
error on K^o induced by the uncertainty on Qoo, the un- 
constrained density dependence of the EoS. Kc and Xc 
are fixed by the existence of a crossing point, and Mc 
is extracted from the correlation analysis based on the 
experimental Eqmr- It is therefore clear that the un- 
certainty on Koo is related to the lack of knowledge on 
the density dependence of the equation of state, repre- 
sented in the present analysis by the skewness parameter 
Qoo- Taking a typical uncertainty for Qo^ of ±400 MeV 
(Table |TT| , Eq. p8)) provides a variation on Koo about 
±40 MeV, compatible with the present uncertainty on 

ifoo. 

In conclusion, the relation Eq. (ITS)) clearly shows that 
the uncertainty on the incompressibility modulus Koo is 
mainly related to that on the quantities Qoo- The reduc- 
tion of the error bar on Ko^ is therefore mostly related 
to a better knowledge of the skewness parameter Qoo , for 
which new experimental constraints shall be found. On 
the contrary the density dependence around the crossing 
point Mc can be more directly constrained from -Egmr, 
as it will be showed below. 
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FIG. 6. Koo as a function of Mc for various Skyrme, Gogny 
and relativistic functionals 



IV. A SIMPLE EXPRESSION RELATING Egmh 

AND Mc 

To provide a complementary view to microscopic ap- 
proaches p]-[ig|. it may be useful to derive an analytic 
relationship between the GMR centroid in nuclei and the 
quantity Mc , in order to enlighten and confirm the results 
obtained with a fully microscopic approach, see Sec. |TT] 
and Ref. [10| . In this section we aim to derive an analyt- 
ical relationship between the centroid of the GMR and 
the relevant quantity of the EoS, Mc Eq. ([7]). 

The energy centroid of the GMR is used to define the 
incompressibility in nuclei Ka :lf|: 



E, 



GMR 



Ih^KA 



(19) 



In order to derive an analytical relationship, (r^) can be 
approximated by 3R'^/5 [323 where R ~ 1.2A^/^ is the 
nuclear radius, yielding: 



Egmr 



h 5Ka 



R V 3771 



(20) 



We shall derive an analytic relation between Ka and 
Mc using the LDA, in order to check, in a complemen- 
tary way to microscopic approaches, the role of Mc in 
determining the centroid of the GMR. 

The following step consists in dividing Ka into a nu- 
clear and a Coulomb contributions, as 

Ka = i^Nuci + i^coui ■ Z'A-4/3 (21) 

where, in the liquid drop approach, i^Nuci is defined as 



K 



Nucl 



^ Koo + K,^,tA-^/^ + Kr 



N - Z 
A 



as in the Bethe Weissacker formula for the binding en- 
ergy [ij . The accuracy of this approach can be enhanced 
with the inclusion of higher order terms [3J| . The quan- 
tities Koo , ^surf and Kr are however poorly constrained 
by the relative small data [i|,[3J]. We prefer instead to ex- 
tract Kj^uci from the LDA which has the advantage that 
i) it was proven to be a good approximation of the mi- 
croscopic calculation |12l | , and ii) the consistency between 
the value obtained for Ka and the Skyrme functional is 
guarantied. 



A. The local density approximation (LDA) 



The nuclear contribution -R'nucI is related to the density 
dependence of the incompressibility K{p) as [13 |. 



K 



Nucl 



A 



d^i 



K{p{r)) 
p{r) 



(22) 



Eq. ([22]) allows to perform the LDA by consider- 
ing the density profile of nuclei, PA{r), in Eq. ([5]), 
where p = pAJf). The LDA give accurate estimation 
of iiTNuci uM- It should be noted that in Eq. (f^ . the 
value of K{p) at saturation density (i.e. Koo) doesn't 
have any specific impact on the Ka value and therefore 
nor on the prediction of E'gmr- Further, due to the ex- 
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h I 207r 
'r I 3mA 



Po/2 



dn 



Po 



- 1 



R 



istence of the crossing area {pc,Kc)— {0.1 fm~^,40 MeV), 
which takes into account both the Skyrme EDFs cross- 
ing (Table m and the relativistic one (Fig[T]), K{p) can 
be approximated to the first order around the crossing 
point by: 



K{p) 



McP 



3 



K, 



(23) 



where M^ is related to the first derivative of the incom- 
pressibility, Eq. (O. 

This first order approximation is relevant as observed 
on the (i?GMR,Afc) correlation of Figs. [2]and[3l Of course 
taking the density dependence of the incompressibility 
as its first derivative around the crossing point remains 
an approximation, which explains the non exactly linear 
{Eqmk^Mc) correlation on Fig. [2] considering Skyrme 
and relativistic EDFs. But still, the correlation among 
the EDF families (Skyrme, Relativistic) are ordered. 

The integral in Eq. ([2^ is taken between po/2 and 
Po, which is adapted to the linear regime around pc and 
corresponds to the typical dispersion of the density val- 
ues around the mean density in nuclei [lO|. Injecting 
expression (|23p in (P^ and assuming a Fermi shape of 
the nuclear density, with diffusivity a ~ 0.5 fm [32|, yield 
the analytical relation between the centroid of the GMR 
and Mc, using Eq. ^ and (pij) : 



fMcP Mc 
\3pc 3 



K)l—^4'ip+'-^Z^A-^A 
/ 1 - p/po p^ 3m J 



1/2 



(24) 



The integral in Eq. ([M]) denotes the nuclear contribu- 
tion whereas the second part comes from the Coulomb 
effects. The Coulomb contribution is evaluated using 
Kcoui=-5.2 MeV d ^. This value is obtained from 
the liquid drop expansion of the incompressibility and 
applied to several Skyrme interactions ^]. It should be 
noted that the Fermi shape is a good approximation of 
the density and we have checked that the diffusivity of 
the density obtained from microscopic Hartree-Fock cal- 
culations (0.47 fm) is very close to 0.5 fm. The use of 
the Fermi density is legitimized by the aim of tracing 
the analytical impact of the quantity Mc on the GMR 
centroid. Equation ((24)) also underlines the important 
role of the quantity M^ on the GMR centroid. On the 
contrary, the incompressibility at saturation density Koo 
doesn't play any specific role in Eq. (IMl) . It is therefore 
rather the quantity M(. which is the relevant quantity to 
be constrained by the GMR measurements. 



B. Results and comparison with the microscopic 
method 



The stability of the results obtained with Eq. (|24|) has 
been studied with respect to the diffusivity value a, the 
LDA prescription (Eq. ([22]) ). the crossing point (pc,Kc) 
values and the integration range. A sound stability is 
obtained against these quantities: the predicted GMR 
centroid doesn't change by more than 10 % by making 
all these variations in relevant physical ranges. 

We first study the behavior of the nuclear contribution 
{.Kcou\—0 in Eq. ([Ml))- Fig. [7] displays the correlation 
between the centroid of the GMR and the Mc value us- 
ing Eq. dMI, for nuclei with A=208 and A=120. A good 
qualitative agreement is obtained with the fully micro- 
scopic results (see Fig. [2]and[3]) in view of the approxima- 
tions performed to derive Eq. (f24|) . The ^-dependence 
is also well described. These results confirm the validity 
of the present approach, and emphasize Mc as a relevant 
EoS quantity to be constrained by the GMR measure- 
ments. It also qualitatively agrees with the microscopic 
results. 

In order to perform a more quantitative study, the 
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with the LDA and including the Coulomb effects (solid lines, 
see text) as a function of Mc- The values for various EDF 
obtained microscopically (squares and dots) are also displayed 
for comparison. 



curves on Fig. [S] displays the centroid of the GMR in 
2°®Pb and ^^°Sn nuclei predicted using the full Eq. ([M)) 
with both the nuclear and the Coulomb contributions. 
The comparison with the microscopic results using vari- 
ous EDFs is also shown. A good agreement is obtained 
in both cases, strengthening again the present analytical 
LDA approach, and emphasising the role of M^. Com- 
paring Figs. [7] and [3 the Coulomb effect on the GMR 
centroid can be evaluated to be about 1 MeV in heavy 
nuclei. 

The almost linear correlation between i?GMR and M^, 



observed on Fig. [8] can be further investigated. Eq. (p4] 
can be rewritten as: 



1/2 



Egmb. = (a(Apo)Me + /3(A^,Po)) ' (25) 



with 
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\P / 
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aKc 



pI 



- P/Pf) P- 
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(27) 



Fixing Pc and Kc, the coefficients a and /3 only depend 
on the nucleus' mass and charge {A, Z), and on the sat- 
uration density po, which is constrained by the charge 
radii. Typical values are a=0.12 MeV and /3=42 MeV^ 
in the case of ^"Spb and a=0.16 MeV and /3=75 MeV^ 
in the case of ^^OSn. 

It should be noted that the LDA approximation (p5)) 
of Egair can be obtained because of the existence of the 
crossing point. In Eq. (|25p . the energy of the GMR de- 
pends on the functional mostly through the parameter 
Mc- In conclusion, the LDA allows to obtain expres- 
sion (|25p relating Eqmr with Mc in a simple and accurate 
form. 

Introducing (Mq^Eq) as the reference point, where 
Mo =1200 MeV and Eq is the corresponding GMR en- 
ergy, one can go one step further and linearise Eq. (j25p 
with respect to Mc — Mq, as 



E^ 



GMR 



2Eo 



-Mr 



En 



aMo 
~2E^ 



(28) 



This is justified for the typical values of Mc, ranging be- 
tween 1000 MeV and 1500 MeV as shown on Fig. [S] 
The almost linear correlation between Eqmr and Mc ob- 
served on Fig. |5]is therefore understood by the present 
approach (Eq. (|28l) ). It clearly shows that the mea- 
surement of the GMR position constrains Mc, which is a 
first information on the density dependence of the incom- 
pressibility. It should be recalled that such a quantitative 
description is not possible with K^o because there is no 
crossing point of the incompressibility at saturation den- 
sity: Eq. ([24l) is not applicable in that case. 



V. CONCLUSIONS 

The relationship between the isoscalar GMR and the 
equation of state raises the question of which EoS quan- 
tity is constrained by GMR centroid measurements. The 



incompressibility modulus Koo alone may not the rel- 
evant one nor the most direct because the more gen- 
eral density dependence of the incompressibility should 
be considered. A crossing area is observed on K{p) at 
Pc — 0.1 fm^"^ among various functionnals. Using a mi- 
croscopic approach, such as constrained-HFB, the slope 
Mc of K{p) at the crossing density can be directly con- 
strained by GMR measurements. This shall assess the 
change of the method in extracting EoS quantities from 
GMR: Mf, is first constrained, and an approximate value 
of Koo can be deduced in a second step [101 • 

The stability of Kc has been demonstrated in the case 
of Skyrme EDFs. A general relationship between Mc and 
Koc is obtained, showing the contribution of the uncer- 
tainty in the density dependence of the EoS which has 
been casted into the quantities Qoo- The Koo value can 
be determined in a second step from the knowledge of 
the Mf, value, requiring a better constraint on the skew- 
ness parameter Qoo, being the main uncertainty for the 
density dependence of the incompressibility between the 
crossing density and the saturation density. One should 
recaU that the Koo value remains 230 ± 40 MeV (17% un- 
certainty) , whereas the quantity Mc is better constrained 



to be Mc=1100 ± 70 MeV (6% uncertainty) [T3|. A bet- 
ter knowledge of higher order density dependent terms of 
E/A{p), e.g. the skewness parameter Qoo, shall help to 
more accurately relate the parameter Mc to the incom- 
pressibility modulus Koo- 

Using the LDA approach and an analytical approxima- 
tion of the density, the microscopic results have been con- 
firmed: the measurement of the centroid of the isoscalar 
giant monopole resonance constrains the first derivative 
Mc of the incompressibility around the crossing point 
Pc — 0.1 fm^-^. A analytical relation between the cen- 
troid of the GMR and the quantity Mc is derived and 
the predicted GMR centroid are found in good agree- 
ment with the microscopic method. 
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